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SUMMARY 


An aK>roximate indicial lift function associated with circulation has been 
developed for tapered, swept wings in incompressible flow. The function is 
derived by representing the wings with a simple vortex system. The results 
from the derived equations compare well with the limited available results from 
more rigorous and complex methods. 

The equations, as derived, are not very convenient for calculating the 
dynamic response of aircraft, parameter extraction, or for determining frequency- 
response curves for wings. Therefore, an expression is developed to convert the 
indicial response function to an exponential form which is more convenient for 
these purposes. The exponential form is nearly as accurate as the form derived 
from the simplified vortex system. 


INTRODUCTION 

The indicial lift function following a unit step increase in wing angle of 
attack is fundamental to the calculation of airplane transient motion. The orig- 
inal derivation of this function (generally designated by k*|(s)) was developed 
in 1925 for two-dimensional incompressible flow by H. Wagner (ref. 1). Since 
that time the indicial response has been calculated for a number of specific 
wing planforms (refs. 2 to 4, for example) and vortex lattice methods are avail- 
able for calculating indicial lift of arbitrary wing planform (ref. 5, for exam- 
ple) . In addition, techniques are currently available for estimating the lift 
on wings undergoing oscillatory motions. (See refs. 6 to 9, for example.) Ref- 
erence 10 contains a review and bibliography on work in unsteady aerodynamics. 

Although available methods appear to be adequate for determining transient 
lift loads, they are generally complex and do not appear to be suitable for 
inclusion in aircraft motion studies (except for a few specific wing planforms) . 

The purpose of this paper is to develop a simplified method for determining 
the indicial lift function associated with circulation. The part associated 
with the "apparent additional mass" is not treated (eq. (31) of ref. 2). A 
second purpose is to provide an elementary approach to unsteady aerodynamics 
which, although it lacks the mathematical elegance of more exact theories, it 
does give the novice a starting point for study in an important and interesting 
^ area of aerodynamics. 


SYMBOLS 


A 


aspect ratio, b^/s 


fa2rbi ,b2 


constants 


b 


ACL(t) 



c 


Cl (0) 

Acj(t) 

Cr 

ct 

i 

K 

k-| (s) 

L 

I 

P,Q,R 

S 

s 

t 

u 

w 

Xo 


YfZ 

a 


wing span 

three-dimensional indicial lift coefficient, a function of time 
steady-state lift curve slope 
wing chord 

section lift coefficient at t = 0 

two-dimensional indicial lift coefficient, a function of time 
wing root chord 
wing tip chord 

= ^P^ 

nondimen sional constant 
Wagner function 
wing lift 

section lift, per unit span 
terms defined in equation (1 2) 
wing area 
Laplace variable 
time 

free-stream velocity 

downwash velocity due to vortices 

distance from three-quarter chord point of root chord to starting 
point of shed vortex along root-chord line 

constants in exponential form of indicial response function 

angle of attack 


dDL 

a = — 

dt 


r total strength of vortex 

Fi vortex strength of ith vortex (i = 1 , 2, . • . , n) 
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A 


sweep of wing quarter-chord line 


As sweep of shed-vortex line 

X wing taper ratio, ^t/^r 

p air density 

T dummy variable 

0 ) angular frequency 


ANALYSIS 

The basic concepts used in this study involve the use of lifting-line 
theory to model the change in lift following a unit step increase in angle of 
attack* The principles are derived for a two-dimensional wing and then extended 
to swept, tapered wings* 


Two-Dimensional Wing 

The two-dimensional wing is represented by a bound vortex located at the 
quarter-chord line* (See fig* 1(a)*) The boundary condition of no flow through 
the wing is satisfied only at the wing three-quarter chord line* If the wing 
angle of attack a is suddenly given a unit increase, the bound (lifting) vor- 
tex strength increases by a value F. Since the total strength (or circulation) 
in the flow field must remain at the value it had before the unit increase in a, 
a counter vortex of strength (-f) is formed near the wing trailing edge* The vor- 
tex strength is such that the normal velocity component at the wing three-quarter 
chord line, induced by the two vortices, is equal and opposite to the component 
of the free-stream velocity normal to the wing surface at the same point* As 
time increases, the shed vortex moves downstream, but retains its original 
strength* This movement would cause a decrease in downwash at the three-quarter 
chord point, but the decrease is offset by the shedding of additional vortices* 
(See fig* 1(b)*) This means, in turn, that the bound lifting-line vortex 
strength increases in order to keep the total circulation unchanged* Conse- 
quently, there is a distribution of vorticity behind the wing* The vortex sheet 
is extending downstream at a rate equal to the free-stream velocity* 

Wagner (ref* 1) obtained a solution for the lift following a unit step 
increase in a by applying a conformal transformation and satisfying the physi- 
cal principle that the velocity at the wing trailing edge is finite at all time 
(Kutta condition)* He derived for the lift, as a function of nondimensional 
distance s = Ut/(c/2), the expression 

2 = 27TpUwk-| (s) (1 ) 


Wagner did not derive an explicit analytical expression for k*| (s) but gave 
only numerical values* Some years later (1936), Kiissner (ref. 11) derived a 
long, slowly convergent series for k*| (s) • 
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The mathematical complexities of obtaining values or expressions for k-j (i) 
arise because of the physical model depicted in figure 1 (b) , or its equivalent 
in the transformed-plane method of reference 1 . The present report develops a 
simpler model and a simple, yet reasonably accurate, analytical expression f6r 
k-) (s) . The shed vortex sheet is replaced by a single shed vortex having the 
same (time varying) strength as the bound vortex and moving downstream at a 
velocity KU. The downwash velocity induced at the three-quarter chord point 
by the two vortices (bound and shed) is easily found by use of the Biot-Savart 
law (or see ref. 12, p. 127), and is 


r /I 

w = + 

2i\\c/2 


Xq + KUt/ 


Solving for T results in 


r = 27TW - 

2 


KUt 

c/2 


1 + 


(2) 


(3) 


The boundary condition is 
w = U sin a 
or, for small angles 
w = ua 

Equations (3) and (4) can 
lift 


met by 


(4) 

be used with the Kutta-Joukowski equation for section 


I = pur 


(5) 


to obtain an expression for the section lift coefficient due to a unit step 
increase in angle of attack 


Ac 2 (t) = 2TT 


/ 


1 - 



5^0 

KUt 

1 + — 

+ / 

c/2 

c/2/ 

Lons (t 

«) , Z 


The distance Xq is related to the starting lift; that 
C2(0). From equation (6), 


(6) 


is the correct value, 
is, Xq is related to 


Ac2 (0) = 2-nil 


\ 


1 \ 



(7) 
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Using Wagner's conclusion that the starting lift is one-half of the steady- 
state lift leads, for the two-dimensional wing, to 


c/2 


= 1 


Equation (6) , therefore, becomes 


Ac 2 (t) = 2tt I 


/ 


1 - 


2 + 


KUt 


c/2/ 


( 8 ) 


(9) 


The term in parentheses is the counterpart of the Wagner function k-| (s) • This 
term is found to approximate the Wagner function very closely if K equals 1/2. 
Equation (9) then becomes 


Ac2(t) 



1 \ 

1 Ut 

2 + 

2 c/2 


( 10 ) 


Garrick (ref. 13) arrived at exactly this same expression by curve-fitting the 
Wagner function. 

The indicial response is approximated closely by use of the single bound- 
vortex line of figure 1 (c) and the shed vortex (which represents the shed distri- 
bution of vorticity) starting at a distance c/2 behind the wing three-quarter 
chord line and moving downstream at one-half the free-stream velocity. 


Three-Dimensional Wing 

The simple vortex model of the two-dimensional wing compares well with 
Wagner's exact analysis; therefore, the same ideas have been tried on three- 
dimensional tapered, swept wings. There are several choices to be made rela- 
tive to the mathematical model, and these choices include: 

(1) At what point should the boundary condition be satisfied? 

(2) Where should the shed vortex be started relative to the wing? 

(3) How should the shed vortex be oriented relative to the wing? 

In order to keep the model as simple as possible, it was decided to satisfy 
the boundary condition only at the three-quarter chord point of the root chord. 
The other two choices are made in a later section of this report. 


5 


The vortex model for the three-dimensional wing is indicated in figure 2^ 
The downwash at the control point can be derived by using the Biot-Savart law 
(or see p. 127 of ref. 12) and is of the form 

r 

w = -(p + Q + r) (11) 

TT ' ' 


where 
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The PrQfR expressions are associated with the bound vortex, the wing-tip 
trailing-vortex pair, and the shed vortex, respectively. The geometric rela- 
tionship b/c^ = A(1 + X)/2 can be used to put equations (12) into more con- 
venient forms 


A (1 + X) 


bP 


A (1 + X) 


2^ ti - tan A 


+ tan A\ 


AH + X) 


sec A - sin A 




+ cos 2 A 


bQ = 


Xo 

Cr/2 2 Cr/2 


1 Ut A (1 + X) 

+ tan A- 


1 - 


A (1 + X) tan* A 


^ 1 ut A (1 + X) 

2 

"a (1 + X)“ 

' J 

r A(i + X) 

1 tan A 

2 

"A (1 + X)" 

+ + tan A~ 

+ 

\ 

+ 


_Cr /2 2 Cr /2 2 


2 

V 

2 


_ 2 


) (13) 


A (1 + X) 


bR = 


A (1 + X) 


sec^ Ag + 


^ 1 Ut \ 

+ j tan Ag 

Cr/2 2 Cj./2y 


1 ut 
7^2 2 ^ 


A (1 + X) 


. / ^ T Ut , 

Ag + + sin Ac 

^Cr/2 2 Or/2 


-i 2 


- tan Ag > 


-o 1 Ut 

— + — cos A <5 

Cr/2 2 cr/2y 


The lift of the wing is determined from the Kutta-Joukowski equation 

L = puFb (14) 

Nondimensionalizing and using equation (11) results in 


CL(t) 


2iTwb 

US(P + Q + R) 


(15) 


The boundary condition requires that w/U = a; hence, equation (15) becomes 


CL(t) 


2irAa 

b(P + Q + R) 


(16) 
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Under steady-state conditions (t , equation (16), with proper substitu- 
tion of equations (13), results in 


2itA 



Solving for 2 tta and substituting the result into equation (16) results in, 
for a unit step in a. 


bP + bQ + bR 


As t becomes very large, equation (18) reduces to 
iCL(t - .) . 

However, equation (17) is not accurate because of the simple vortex-system rep- 
resentation of the wing. A concept proposed is to use the best available source 
for (^Lct) rather than equation (17). In this manner, attaining correct val- 
ues for equation (18) is insured, at least for large values of time. 

Before equation (18) can be used (with eqs. (13)), the terms Xo/(Cr/2) 
must be determined, and a choice made for the sweep angle of the shed vortex. 

As noted previously, for a two-dimensional wing, a value of Xo/(Cr/2) =1.0 
is appropriate. This value is retained for simplicity so that the inboard ends 
of the shed vortex start at a distance Xq/(Ci-/2) = 1 behind the three-quarter 
chord point of the root chord. 

Selection of the sweep of the shed vortex was made by trying reasonable 
values; that is, (1) parallel to the wing trailing edge, and (2) parallel to 
the wing quarter-chord line. Comparison with results from more exact methods 
showed that the shed vortex sweep angle equal to that of the wing quarter-chord 
line was the better choice. 
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Based on these choices, equations (13) become 


A(1 + X) 


bP 


A(1 + X) 


^2 A - tan A 


2 1 I 

fA(l + X) “I 

2 

\ 

sec A sin A 

+ cos 2 A 


2 



+ tan A; 


bQ 


A(1 + X) 1 Ut 

1 + tan A + 

2 2 Cj./2 


AO + X) 

1 - tan A 


1 ut A(1 + X) 

1 + _ + tan A 

2 Cr/2 2 


A(1 + X) 


A(1 + 

1 

2 


X) 


tan A 


' 2 

1 

"ao + X) 

-r 

2 


/ (19) 


A(1 + X) 


bR = 


A(1 + X) / 1 Ut \ 

sec2 A + 1 + - tan A 

2 \ 2 Cr/2/ 


1 Ut 
^ + - 


2 Cr/2 


A(1 + X) 


/ 1 ut \ 1 

2 1 

r/ 1 ut \ 1 

2 

A + 1 + 1 sin A 

+ 

h + cos A 


V 2 Cc/2; 

i 

2 Cr/2/ 



- tan A 


Equations (18) and (19), with suitable selection of an expression for 



(use of ref. 14 or 15, for example), are the basic equations for the indicial 
response for tapered, swept wings in incompressible flow. Although the equa- 
tions are somewhat lengthy, several factors appear repeatedly and are constant 
during the lift buildup. 


LIFT VARIATION FOR ARBITRARY CHANGE IN ANGLE OF ATTACK 

The lift of a wing undergoing an arbitrary change in angle of attack can be 
determined by use of the indicial response function in Duhamel's integral; so 
that 


CL(t) 



(t - T) 6t(T)] 


dT 


( 20 ) 


Frequency response characteristics can be obtained easily by using the Laplace 
transform of equation (20) so that 

Cl(s) = ACi,(s) d(s) = [ACl(s)] (s)[a(s)] (21) 
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from which 


C£,(iw) = (i(D) [a(iw)J 


Because of the manner in which time enters into equations (18) and (19)/ closed- 
form Laplace transforms are not known to exist for the general indicial lift 
function ^^(t). Therefore, a form for which closed-form Laplace transforms 
do exist, which would also be convenient, has been sought. In reference 2 the 
indicial response for several specific elliptic wings were given in exponential 
series of the form 


ACr. = 


L = ao + ai 




In the present study, equation (18) was approximated by just one exponential 
term, namely. 


^L(t) 




that leads to 


1 

(=I.C.)„ ■ bP + bQ + bR 


bP + 1 + 


. 1 - 


A(1 + X) 

1 tan A 

2 


12 r 12 

A(1 + X) A(1 + X) 

1 tan A + 

2 2 


The parameter y was determined by evaluating equation (24) at 
Ut/(Cr/2) = 0 and solving for y. The resulting equation is 


y = 1 - = (bP + 1 + 

[bP + bQ + bR]ut/(Cr/2)=0 I 


A(1 + X) 

1 tan A 

2 

12 r 12 

A(1 + X) A(1 + X) 

tan A + 

2 2 
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where 


(bP) = 


A(1 + X) 


C A(1 + X) 

sec^ A - tan A 

2 

1 + tan A 


A(1 + X) 


sec A - sin A 




+ cos^ A 


J 


A(1 + X) 

1 + tan A 


(bQ) ut/ (Cr/2) =0 - 


A(1 + X) 

1 + tan A 


A(1 + X) 


A(1 + X) 

1 - tan A 


If A(1 + X) 1 

2 

+ 

< 

1 

2 

/ 1 - tan A 

+ 



L ^ J 


2 



A(1 + X) 


A(1 + X) 


sec 


2 A + tan A 


(bR)ut/(Cr/2)=0 = ; 




=7 

> 

+ 

1 

2 

sec A + sin A 

2 

+ cos^ A 


- tan A 






> (26) 


J 


Determination of the parameter z is more difficult. If equation (23) is 
solved for z, the result is 


1 


z = - 


Ut/ (Cr/2) 


r 

-Ti 

1 

AC^Ct) 


lo9e - 

1 - 




( 

L 


J 


(27) 


n 



Evaluation of equation (27) at Ut/(Cj./2) = 0 yields the indeterminant form 

f ACl(O) "I 

0/0, since y is precisely equal to 1 - . However L 'Hospital's rule 

(%)ss_ 

can be used with equation (27) and the resulting expression evaluated for z. 

By using L 'Hospital's rule 




X 



bP + bQ + bR 

-Ut/(Cr/2)=0 
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or 


z = 


1 , 

- -/bP + 1 

y ' 




A(1 + X) 

1 - tan A 


1 - 


A(1 + X) 


tan A 


A(1 + X) 


5> 




(bP + bQ + bR)2 


Ut/ (Cr/2) =0 


Ut 

<^/2 


(bP + bQ + bR) 


(28) 


Ut/(Cr/2)=0 


The derivatives appearing in equation (28) were evaluated, and the results are 


d(bP) 


'(^i) 


bt/(Cr/2)=0 


d(bQ) 




d(—\ /T, .1^ r*<’ + 

11' * "" *J • L“^~J J 


d<bR) 

1 A 

“ - ~<‘»R)ot/{Cr/2)=0 + - 

(1 + X) tan A 

Ka^)J 

IfAd + X) o 

,/ sec A + sin A + cos^ A 

Ot/(Cr/2)-0 

VL 2 J 


A(1 + X) 




X) 


sec^ A + tan A 
A(1 + X) 


A(1 + X) 


sec A + sin 


aJ sin A + cos2 ^ 


{ 


-12 -^3/2 

sec A + sin A + cos^ A> 


(29) 
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Equations (25) and (26) were used to calculate y for a wide range of wing 
sweep angles, aspect ratio, and taper ratio. Equations (28), (26), and (29) 
were used to calculate z for the same range of geometric variables. Results 
from the calculations are presented in figure 3 and can be used in equation (23) 
for calculating indicial response. 


COMPARIS(»J OF RESULTS WITH OTHER METHODS 

This section compares results from the methods developed in this study with 
other methods for which numerical results are published. 


Two-Dimensional Wing 

The two-dimensional wing is a rather special case having A ^ X = 1 , 

A = 0. The indicial-lift function given by the equation is relatively simple 
(eq. (10)) and does possess a closed-form Laplace transform. There is no need, 
therefore, to use the exponential form for the indicial lift. Results from 
equation (10) are compared with those of references 1, 2, 4, 5, and 11 in fig- 
ure 4. The results from equation (10) are in excellent agreement with those of 
references 1, 2, 4, and 11, which are exact solutions or numerical approxima- 
tions of exact solutions. The success of equation (10) in matching these exact 
solutions is attributable to the selection of Xq (the starting point of the 
shed vortex) and the velocity of the shed vortex to match the Wagner solution. 


Three-Dimensional Wings 

Vortex systems and doublet distributions used to model three-dimensional 
wings are generally rather complex and lead to involved calculations for deter- 
mining indicial lift. Numerical results, therefore, are available for only a 
few specific wings for comparison with equation (23) . 

Rectangular wings .- Rectangular wings have X = 1.0 and A = 0. Results 
from use of equation (23) , with values of y and z from figure (3) , are com- 
pared with results from other more exact (and more complex) methods in figure 5 
for wings of aspect ratio 6 and 4. The figure shows that the approximate method 
of the present paper is in good agreement with the results of references 4 and 5. 

Triangular wings .- Triangular wings have X = 0 and, in addition, the wing 
geometry leads to a simple relationship between aspect ratio and the angle of 
the quarter-chord line, that is, 

3 

tan A = - (30) 

A 

Results from use of equation (23) and values of y and z fran figure 3 are 
compared with results of reference 16 in figure 6 for triangular wings of aspect 
ratio 4 and 2. In these cases the agreement between the results of the present 
approximate theory with results of reference 16 is excellent. 
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Limiting case with A ^ The parameter y for wings having aspect ratio 
of zero can be determined readily frcan equations (25) and (26) and is 


= 0 

The indicial response function, therefore, is 
’ACL(t) 




ss 


= 1.0 


A>0 


This result agrees with the results of reference 16. 


(31) 


(32) 


CONCLUDING REMARKS 

An approximate indicial lift function associated with circulation has been 
developed for tapered, swept wings in incompressible flow. The function is 
derived based on representing the wings by a simple vortex system. Cmparison 
of results from the derived equations compare well with the limited available 
results from more rigorous and ccxnplex methods. 

The equations, as derived, are not very convenient for calculating the 
dynamic response of aircraft, parameter extraction, or for determining frequency- 
response curves for wings. An expression, therefore, is developed to convert 
the indicial response function to an exponential form which is more convenient 
for these purposes. The exponential form is nearly as accurate as the form 
derived from the simplified vortex system. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
May 5, 1978 
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(b) Vortices some finite time after a step increase in a. 
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(c) Simplified vortex system used in this study. 

Figure 1 .- Vortex representation of a two-dimensional wing following a step 

increase in angle of attack. 
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Figure 5.- Concluded. 
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(a) Aspect ratio 4.0. 

Figure 6.- Indioial lift for delta (triangular) wings 
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